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ABSTRACT 

In  this  paper,  shape  change  of  the  elastic  and  rigid  boundary  in  a  micro-channel  is  simulated.  The  rings  are 
considered  as  elastic  boundaries  immersed  in  the  fluid  flow.  In  this  study,  we  use  the  Immersed  Interface  Method  to 
simulate  the  flow  around  rigid  objects  and  movement  of  the  flexible  objects.  In  this  method,  there  is  no  compulsion  to 
match  the  object  mesh  and  the  flow  mesh.  The  results  were  found  to  be  in  good  agreement  with  available  data.  First,  the 
simultaneous  motion  of  multiple  membranes  in  a  micro-channel  and  their  interaction  with  each  other  and  with  flow  are 
evaluated.  In  addition,  simultaneous  presence  of  multiple  membranes  would  result  in  a  reduction  in  the  flow  velocity. 
Finally,  the  effects  of  changes  in  the  elastics  and  bending  moduli  on  the  deformation  of  cell  are  considered. 
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INTRODUCTION 

The  immersed  interface  method  (IIM)  is  the  most  appropriate  approach  for  solving  the  problem  of 
fluid-solid  interaction.  The  basis  of  IIM  is  to  add  a  forcing  term  as  a  source  term  to  the  Navier-Stokes  [1].  Feng  and 
Michaelides  [2]  were  the  first  to  IIM  and  simulated  suspensions  of  rigid  disks  in  2D.  Le  and  Zhang  [3]  used  in  their 
work  a  hybrid  LBM-  IIM  and  noticed  that  the  computed  velocity  profiles  can  deviate  greatly  from  theoretical  ones 
even  for  very  simple  flow  situations,  both  in  the  immersed  boundary  layer  and  the  bulk  region.  Dupuis  et  al.  [4] 
studied  how  the  coupling  method  of  the  forcing  term  between  the  Eulerian  and  Lagrangian  grids  could  affect  the 
results  for  the  flow  over  an  impulsively  started  cylinder  at  moderate  Reynolds  (Re)  number.  Wu  and  Shu  [5] 
proposed  a  new  version  of  IIM,  which  could  well  consider  the  effect  of  external  force  on  the  momentum  flux  as 
well  as  the  discrete  lattice  effect.  JiSeok  and  SangHwan  [6]  presented  a  numerical  scheme  for  fluid-structure 
interaction,  especially  for  elastic  structures.  They  employed  a  hybrid  IIM  using  an  improved  direct  forcing  scheme 
for  the  fluid,  and  a  finite  element  method  with  Euler  beam  elements  for  the  elastic  plate.  Zhang  et  al.  [7,  8]  also 
used  a  combination  of  the  IIM  and  the  LBM  to  investigate  the  microscopic  hemodynamic  and  hemorheological 
behaviors  of  discrete  RBCs  in  shear  flow.  They  noted  that  three-dimensional  simulation  of  RBCs  is  required  to 
attain  accurate  results.  In  the  present  work,  first,  the  simultaneous  motion  of  multiple  membranes  in  a  micro- 
channel  and  their  interaction  with  each  other  and  with  flow  are  evaluated.  Finally,  the  effects  of  changes  in  the 
elastics  and  bending  moduli  on  the  deformation  of  cell  are  considered. 
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GOVERNING  EQUATIONS 

Navier-Stokes  equations  for  the  incompressible  flow  are  as  follows: 

F‘U  =  0  (1) 

p  +  u  •  =  —  Vp  +  fjF2  u  +  f1  (2) 

With  boundary  and  initial  condition  [9]: 

«Ue  =  uh  (3) 

aU,  0)  = 

Where,  u  is  the  fluid  velocity,  p  is  the  pressure,  p  is  the  density,  and  p  the  viscosity  of  the  fluid.  The  force  F  which 
has  the  form 

F&  i )  ='  t)5(x-$(s,  t)}ds  (4) 

The  motion  of  the  interfaces  is  shown  as  below  formulation 

u{x(s,t),  t)  =  =  j£  u(x,t}5  -  ,YC>,  t)}  dx  (5) 

The  force  density  using  an  expression  of  the  form, 

/(s,  t)  =  K  (^ts)  -  £ Cs,  o)  g 

The  force  is  expressed  as 

/  (s,  t)  =  a e/(s,  tr)  dt r  -(-  fl  $(s,  t)  (7) 

The  respective  normal  and  tangential  components  of  the  force  density  /  =  f(s,  t'},  /7  and  f,  =  f(s,  t').  t  can  be 
related  to  the  jump  conditions  forpressure  and  velocity  as  follows  [8]: 

[m]  =  0,  =  —f2i,  1  =  0  (8) 

W=A  .W=S  (9) 

=  k/t,  n  (10) 

In3ff]  =  -[nitpj  t  b?f]n  t  t  ,i?[wf]w.it  (11) 

We  consider  the  below  relations  as  a  simple 

[“*]  =  [“J«l+  [“jpki  (12) 

[ity.y]  =  [%]«s  +  2[i*fJ^T!  t  [^]t|  (13) 
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RESULTS  AND  DISCUSSIONS 

Validation 

In  this  section,  we  investigate  the  interaction  of  the  flow  and  a  circular  ridge  cylinder  located  at  the  center  of  the 

channel.  Vortex  shedding  frequency  is  /v.  The  Strouhal  number  (St= - )  and  the  average  drag  coefficient  are  compared 

with  other  results  in  Table  1,  which  is  in  good  agreement.  It  can  be  seen  that  with  the  increase  in  Reynolds  number,  the 
average  drag  coefficient  decreases,  but  the  Strouhal  number  increases.  Increasing  the  Reynolds  number  increases  the  rate 
of  the  vortices  spreading.  That  is,  in  a  certain  time,  more  vortices  are  spread  in  the  flow. 

Table  1  Comparison  of  the  present  results  with  other  results 


Cp 

St 

Braza  et  al.  [10] 

1.359 

0.16 

Re=100 

Ji  et  al.  [11] 

1.376 

0.169 

present  result 

1.380 

0.16 

Braza  et  al.[10] 

1.386 

0.20 

Re=200 

Ji  et  al.  [11] 

1.354 

0.20 

present  result 

1.355 

0.20 

Simultaneous  Motion  of  Three  Membranes 

In  this  section,  effect  of  membrane’s  location  on  each  other  is  investigated,  Three  membranes  were  placed  on  the 
micro-  channel’s  axis  of  symmetry  (Figure  la)  in  Figure  lb,  deformation  of  the  membrane  at  time  t  =  0.5  is  observed. 
It  is  observed  that  pressure  and  velocity  at  top  and  bottom  of  membrane's  center  had  symmetry.  So,  there  are  no  pressure 
differences  in  the  vertical  direction  so  as  to  deviate  membrane  from  the  center  of  the  micro-channel.  The  low  pressure 
difference  created  around  the  membranes  at  x-direction  would  create  a  force  from  fluid  to  membrane  in  fluid  motion 
direction  from  left  to  right.  As  can  be  seen,  all  three  membranes  would  gain  more  convexity  in  front  and  more  concavity 
from  behind,  so  as  to  overcome  the  flow  resistance.  This  result  can  also  be  seen  in  experimental  results  of  [12,  13]. 


3  2  l 

ooo 


(a) 


DDD 

i - i - i - 1 


(b) 

Figure  1:  (a)  Initial  Position  of  the  Three  Membranes, 
(b)  Their  Motion  and  Deformation  at  Time  T  =  0.5 
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The  Effect  of  Increasing  the  Reynolds  Number  on  the  Flow  Pattern 

The  problem  geometry  in  this  case  is  shown  in  Figure  2.  All  dimensions  of  the  cylinders  are  dimensionless  by 
the  length  of  the  square  obstacle.  The  length  of  the  flow  input  (the  distance  between  the  front  edge  of  the  square 
cylinder  and  the  channel  entrance)  is  L  /  5.  The  Reynolds  number  is  Re  <300.  Figure  2  shows  the  flow  lines  for  several 
different  Reynolds  numbers.  When  Re  =  1,  due  to  the  predominance  of  Viscose  forces,  a  steady  flow  will  be  formed 
without  any  separation.  As  the  Reynolds  number  increases,  separation  of  the  flow  occurs  in  the  back  of  the  barrier. 
The  onset  of  the  separation  phenomenon  and  the  formation  of  stationary  vortices  for  the  pre-sent  problem  is  Re  =  3 
(Figure  2a).  As  can  be  seen,  this  value  is  for  a  square  barrier  smaller  than  the  corresponding  value  for  the  flow  around 
the  circular  cylinder.  The  reason  for  the  smaller  value  of  this  value  in  the  flow  around  the  square  barrier  can  be  related 
to  the  existence  of  sharp  edges  in  this  problem.  When  the  Reynolds  number  of  the  upstream  flow  reaches  the  critical 
Reynolds  number,  alternating  unsteady  vortex  phenomenon  occurs.  The  critical  Reynolds  number  for  the  studied  issue 
is  calculated  with  blockage  ratio  of  0.125,  about  55.  References  [14,  15]  have  provided  54  and  70  respectively  for 
critical  values  of  Reynolds.  Figure  4b  shows  the  stream  lines  for  Re  =  200.  Since  the  flow  in  this  case  is  unsteady,  the 
shape  of  the  stream  lines  is  drawn  in  a  special  moment.  Flow  in  Re  =  200  fn  addition  to  the  vor-tex  flow  behind  the 
cylinder,  smaller  vortex  flows  will  be  formed  adjacent  to  the  lateral  edges  of  the  cylinder.  fn  other  words,  the  separation 
region  from  the  back  of  the  cylinder  penetrates  to  the  front  edge.  Based  on  the  results  of  the  calculations,  the  penetration 
phenomenon  of  the  separation  region  to  the  front  edge  can  be  seen  in  the  Reynolds  range  greater  than  130. 
This  phenomenon  has  significant  effects  on  the  flow  characteristics  that  will  be  addressed  in  the  following  sections. 


Figure  2:  The  Stream  Lines  for  a)  Re  =  3,  b)  Re  =  200 

Figure  3  shows  the  variation  in  vortex  region  length  (steady  flow  with  symmetric  vortices)  versus  Reynolds 
number.  The  calculations  confirm  the  linear  increase  in  vortex  region  length  in  terms  of  Reynolds  number  increase. 
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Figure  3:  The  Vortex  Region  Length  Versus  the 
Reynolds  Number  in  the  Permanent  Flow 

CONCLUSIONS 

In  this  study,  we  use  the  Immersed  Interface  Method  to  simulate  the  flow  around  rigid  objects  and  movement  of 
the  flexible  objects.  In  this  method,  there  is  no  compulsion  to  match  the  object  mesh  and  the  flow  mesh.  The  results  were 
found  to  be  in  good  agreement  with  available  data.  First,  the  simultaneous  motion  of  multiple  membranes  in  a 
micro-channel  and  their  interaction  with  each  other  and  with  flow  are  evaluated.  The  numerical  simulation  of  the  flow  of 
fluid  through  one  or  a  set  of  objects  that  causes  the  flow  to  separate  from  the  surface  of  them  has  been  the  subject  of 
interest  by  researchers  over  the  past  few  decades.  One  of  the  most  important  types  of  these  objects  is  those  with  a  square 
cross  section,  which  have  important  and  diverse  applications  in  different  industries.  The  flow  around  a  square  cylinder  was 
simulated  and  we  surveyed  different  flow  patterns.  The  changes  in  the  vortex  region  length  and  the  Drag  coefficient  were 
investigated  in  different  Reynolds. 
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